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Abstract 



We show that in the neighborhood of the tripling bifurcation of 
a periodic orbit of a Hamiltonian flow or of a fixed point of an area 
. preserving map, there is generically a bifurcation that creates a "twist- 

^ | less" torus. At this bifurcation, the twist, which is the derivative of 

' the rotation number with respect to the action, vanishes. The twist- 

less torus moves outward after it is created, and eventually collides 
with the saddle-center bifurcation that creates the period three or- 
bits. The existence of the twistless bifurcation is responsible for the 
^ ■ breakdown of the nondegeneracy condition required in the proof of the 

KAM theorem for flows or the Moser twist theorem for maps. When 
the twistless torus has a rational rotation number, there are typically 
reconnection bifurcations of periodic orbits with that rotation number. 
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1 Introduction 



The dynamics in the neighborhood of an elliptic periodic orbit of a two-degree 
of freedom Hamiltonian flow, or equivalently, an elliptic fixed point of an area 
preserving map, can be elucidated by consideration of their formal normal 
forms. When the rotation number, u, of the elliptic orbit is irrational, the 
normal form is called the Birkhoff normal form. Let J denote the transverse 
action (the "symplectic radius") and 6 be its conjugate angle. The Birkhoff 
normal form for the flow is 

H(J)=UJJ+ 1 -T {) J 2 + 1 -T 1 J3 + ... . (1) 

2 o 

For an area-preserving map with an elliptic fixed point, the Birkhoff normal 
form is 

J ^ J 

6 ^ 6 + 2ntt(J) (2) 
where the rotation number is 

Q(J) = uj + r J + ~ti J 2 + . . . . 

This map is also the time 2tc map of the Hamiltonian flow. We define the 
"twist" to be the derivative of the transverse rotation number with respect 
to the action: 

/ T N dfl T 
r(J) = — = T Q + T l J + .... (3) 

In general, the Birkhoff series do not converge. The actual dynamics of 
the system typically includes a chain of n islands located near each radius 
where Q = m/n, a rational number. This chain is constructed from an elliptic 
and hyperbolic pair of periodic orbits of period n. These orbits are born in an 
rz-tupling bifurcation that occurs when uj passes through the rational value 
m/n. 

The character of the bifurcation as uj passes through m/n when tq ^ 
depends upon n [ |dA90| , |MH92| , [Sch98|1 . When n > 4, the elliptic orbit remains 



stable, and the bifurcation gives rise to a chain of n islands. As we recall 
in the Appendix, for n = 4 the origin can become unstable providing the 
resonant term is large enough. When n = 3, the resonant term generically 



2 



dominates the twist term near the origin, and the origin is unstable at the 
bifurcation point. 

When the rotation number Q is a monotone function of J, or equivalently, 
the twist does not vanish, the map is called a "monotone twist map." For 
this case, Aubry-Mather theory applies |Mei92|l and implies in particular that 
there is a pair of rotational m/n periodic orbits for each rational in the range 
of Q. The nonvanishing of the twist also corresponds to the isoenergetic 
nondegeneracy condition required for KAM theory ||Arn78 1 , in the proof of 
the existence of invariant tori. 

When the twist vanishes, the dynamics can be much more complicated. 
For example, "reconnection bifurcations," occur near an extremal point of 
Q [ |HH84| , |HH95| |, and the renormalization operator for the destruction of 
invariant circles gives a distinct universality class for circles that cross the 
zero twist line as opposed to those which do not | |dCNGM96[ |. Recently 
||Sim99|| it has been shown that the "meandering curves" that appear when 
a twistless curve passes through a rational rotation number are stable under 
small perturbations. An extension of standard KAM theory |pdlL98|| shows 
that in two parameter families of area preserving maps a diophantine twistless 
curve persists. 

While it appears that To should generically be nonzero, and so the twist 
is nonvanishing at least in some neighborhood of the origin, we will show 
in this paper that this is not true whenever the rotation number u passes 



through 1/3. This answers the question about genericity raised in [[Sim99 



To show that the twist generically vanishes, we begin with the normal 
form in the neighborhood of a tripling bifurcation. Generally, when the rota- 
tion number of the elliptic periodic orbit is rational, there are resonant terms 
that cannot be transformed away. If we keep only the first such resonant 
term, the Hamiltonian becomes (see e.g. | dA90| . 1\1H92 |) 



H(I, 9,t) = uI + AI 2 + --- + BT n/2 cos(n9 - mt) + . . . . 

We now use the variable / to denote the "action" for this system; it is not a 
true action variable. Here the system is at resonance when u = m/n. The 
time dependence in H can be eliminated by a transformation to rotating 
coordinates, defining = 9 — mt/n. This gives the new Hamiltonian H = 
H — ml /n: 



H(I, 



eI + AI 2 + --- + BI n/2 cos(n0) + . . . 



(4) 
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where e = uo — m/n measures the frequency difference from the resonant case. 
The time 2ir map of H is an area preserving map which has if as a conserved 
quantity. 

When the resonant coefficient B vanishes, then the coordinate I is the 
true action, J, and the twist for Eq. @ is r = 2 A + O(J), as we can see by 
comparison with Eq. (Q). However, when B ^ 0, the action is modified and 
so is the twist. Moreover when n = 3, the resonant term is of lower order 
than the first twist term, and even the 0(J°) terms in the twist must be 
corrected. We will compute the twist for this system in Sec. |2|. 

We also consider the resonant normal form for an area-preserving map in 
cartesian coordinates {p,q)- This is most easily written in complex coordi- 
nates, 



'2Ie 2me =p + iq, (5) 

and in terms of the multiplier of the fixed point 

A = e 2muj . (6) 

Then the normal form can be written 

z' = X(z + iaz 2 z + ■■■ + Pz 71 - 1 + ...). (7) 

where the omitted terms, as we will see below, include terms that are required 
for the map to be area-preserving. By comparison with Eq. (^), the twist, 
when the resonant terms vanish is r = a /it + O(J). The resonant terms, 
however, will modify the twist; indeed when n = 3, the resonant terms are 
of lower order than the first twist term, and so even the 0(J°) terms in r 
should be changed. 

We will show that whenever the resonant term is nonzero, the twist van- 
ishes in the neighborhood of the tripling bifurcation. We will do this by 
assuming that e = u — m/n^0in Eq. (f|) and Eq. (|7p, so that they can be 
transformed to Birkhoff normal form. This gives an expression for the twist 
that diverges at resonance. The twist, however, is well defined away from 
resonance, and we will see that r = at some uo = a;o in the neighborhood 
of 1/3. At this rotation number a twistless torus is created at the origin. 
As lu moves away from u>o, the twistless torus grows, corresponding to an 
extremum in Q(J). 

As Moser showed ||Mos94j| , the twist for the Henon map vanishes at ujq ~ 



0.29. We show in addition that the twistless circle moves away from the 
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origin as u moves towards 1/3. This implies, for example, that when the 
rotation number of the twistless circle passes through a low order rational 
number, such as 3/10, a reconnection bifurcation ||HH84| , [HH95|| should occur. 
The dynamical consequences of this have already been observed in [|W90| , 
[5T981 , though without explaining their genericity. Using our calculation of 
the higher order twist, we will obtain a good approximation for the position 
of this bifurcation. 

More generally, the twist can vanish at any rotation number. For example, 
for a cubic map, we will show that with the choice of two parameters we can 
make To vanish at u = 1/5 or 2/5, which results in the instability of the 
elliptic point at the bifurcation. Similarly the seventh order resonance can 
be generically destabilized in three parameter family of quartic maps. The 
instability of these resonant twistless maps is proved in the appendix. 



2 Resonant Hamiltonian flows 

In this section we will take advantage of the fact that the resonant normal 
form Eq. is integrable, to obtain exact expressions for the twist near 
resonance. We begin by rewriting Eq. (f|) in the form 

H(I,(f>) = P(I) + Q(I)cos(n<p), with 

P(I) = eI + AI 2 + ... (8) 
Q{I) = I n ' 2 (B + ...), n>2, 

where the dots denote finite polynomials in /. Our goal is to transform I to 
the true action variable, J = ^- § I dtp for Eq. @ and obtain the expression 
for the rotation number Q(J). To do this, we assume explicitly that e ^ 
so that P has a first order zero at the origin. 

To find the period, T(h), we solve H(I, <fi) — h for <fi and substitute this 
into the differential equation for / to obtain 



T(h) 



dl 



v/q 2 -(h- py 



Here we have taken into account that the extrema in I are visited n times 
for a full turn of <fi; therefore the integral has been multiplied by n. 

The action is given by the area under the curve I(h,4>), but it is better 
to obtain it by integrating —0(1, h)dl, because otherwise we would need to 
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solve a quartic or higher equation for /. The area under I is given by the 
difference between the maximum area 2irl max and 2n times the area under 
the curve between extrema, since the cosine is even, and each extremal is 
visited n times in one loop. Hence the action turns out to be 



J{h) = — 27r/ max (/i) -2n 0(7, h)dl . (9) 



As usual one can verify that T(h) = 27rJ'(/i), since the boundary terms 
cancel. 

As expected, the period is an Abelian integral of the first kind because 
P(I) and Q{I) 2 are polynomials in I. The action can also be turned into an 
Abelian integral by partial integration of Eq. Again the boundary term 
cancels and we obtain 

Lemma 1. The action J{h) and the period T(h) of the resonant normal 
form Hamiltonian Eq. are Abelian integrals on the Riemann surface 

T:y 2 = R{z) = Q{zf-{h-P{z)f. 

They are given by 

= fP- W9W+ QW-P W )^ 
/ 7 Q{z)y 

T(h) = (f-dz, (10) 
J*y y 

where 7 corresponds to the real cycles ofT. 

The degree of R(z) is at least n because the lowest order term in Q(z) is 
z n / 2 . For n = 3 the integrals are elliptic, for larger n they are hyperelliptic. 
The equilibrium points of Eq. (Q) correspond to the double roots of R(z). 
The values of h for which this occurs are called critical values and they can 
be determined by the discriminant of R(z). The simplest critical value is 
h = 0, since there is always a double root at z = 0: 

R(z)\h=0 = Q(Z) 2 - P(Z) 2 = -Z 2 6 2 (l + ...). 

This corresponds to the fixed point at the origin. Since R"(0) < at h = 
there is a maximum between the colliding roots so this corresponds to a 
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vanishing real cycle 7.Q All that is left from the Abelian integral in this case 
is the residue of the pole at z = 0. The limiting period is therefore given by 

1 2tt 

T(0) = 2vnRes- = — 

z=o y e 

This is a trivial result, because e was designed to be the deviation in rotation 
number from m/n in the first place. Note, however, that the original variable 
/ was not the true action. It is the fact that the resonant term is of order 
n/2 that the frequency at h — stays the same. A similar calculation will 
lead to the twist of the origin, which can be changed by the resonant term if 
n = 3. 

2.1 Tripling bifurcation for flows 

For the tripling bifurcation the Hamiltonian from Eq. (|J) through quadratic 
order is 

H{I,<j>) = el + AI 2 + BI 3/2 cos(3(f)). (11) 

Assuming that both A and B are nonzero, we can eliminate both parameters 
in H upon defining 

A 2 ~ A 3 A 
/ = ^2 7 ' h = J^i h ' £ = J^2 £ > and0 = + V, 

where ip — if AB > and ip = n/3 otherwise. This scaling leaves only the 
two essential parameters (e, h), so that the energy equation becomes 

h = il + l 2 + P /2 cos(30) . (12) 

The corresponding elliptic curve is 

R( z ) = z 3 -(h-£ Z - z 2 ) 2 . 

We consider e as the bifurcation parameter and the energy h is the parameter 
selecting a particular torus. 

The phase portrait in Cartesian coordinates, H((p 2 + g 2 )/2, arctanp/g), 
always has one equilibrium point at the origin corresponding to the line 

^'"(O) ^ 0: for n = 3 it vanishes at e = 1/2, but then R^(0) - 12 < 0. 
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Figure 1: Bifurcation diagram for Eq. (|T^). The horizontal axis corresponds to 
the fixed point at the origin. The upper of the two remaining curves represents 
the saddle 3- fold orbit and the lower one the elliptic 3-fold orbit. At e = 9/32, 
where the two lines meet, these orbits are collide in a saddle-center bifurcation. 
The shaded area is energetically forbidden 



h = in Fig. |1[ The other equilibria are easily found from the zeros of the 
discriminant D of R(z), which is given by 

D(h,e) = -h 3 256/? + (27- 144^+ 128e 2 )h + 4i 3 (-l + As) 

Apart from the triple zero at h — 0, D has two distinct real roots for h when 
e < 9/32. These are the energy levels for one of the 3-fold equilibria which 
correspond to the period three orbits in the original, nonrotating frame. The 
point (e, h) = (9/32, (3/16) 3 ) corresponds to the cusp in Fig. [I], where the 
period three orbits collide. There are four regions in the bifurcation diagram; 
the corresponding phase portraits are shown in Fig. |^. In region (i) there 
is one interval of positive R corresponding to one interval of real momenta, 
hence one torus. Crossing the critical lines to regions (ii) and (iv) creates a 
second interval of positive R out of a double root. Entering region (iii) from 
these regions by crossing another line of critical points destroys a positive 
interval in a double root so again there is only one torus corresponding to 
each point in region (iii). Regions (i) and (iii) can also be left by destroying 
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the one positive interval in a double root so that no motion at all is possible 
(the shaded region). 

The details of which torus in the phase portrait belongs to what region are 
shown in Fig. ^|. Consider the one positive interval of region (i): For e < it 
corresponds to tori sufficiently far away from the origin (part of dotted line); 
for 0<£<l/4it corresponds to all tori outside the separatrix (dotted line); 
for 1/4 < e < 9/32 it corresponds to all tori outside the separatrix (dotted 
line) and also to tori sufficiently close to the origin (part of solid line); for 
i > 9/32 it corresponds to all tori (dotted line). For the other regions the 
correspondance is analogous. The boundaries of the regions, where two roots 
of R coalesce, correspond to the equilibrium solutions. Note that the origin 
is a minimum of H((q 2 +p 2 )/2, arctanp/g) for positive e and a maximum for 
negative e, showing that it is stable when e ^ 0. 




Figure 2: Bifurcation diagram for the Tripling Bifurcation in (i,h). The thick 
curves represent the energy levels of the equilibrium points. The phase portraits 
correspond to fixed values of e, and the thin, vertical lines represent the energy 
ranges for the topologically distinct families of tori, shown in each phase portrait. 
Each group of thin lines corresponds to a single e, but they are drawn slightly 
skewed for clarity. 
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2.2 The twist 



The twist r(h) of a torus with energy h is given by Eq. (|3|), thus 

or, expanding near the origin in a series in J, we have 

r(J) = -e 3 J"(0) + e 4 [3eJ"(0) 2 - J'"(0)] J + 0(J 2 ) 

To compute these first two terms in the twist, we must find the second and 
third derivatives of J(h) at h = 0. 

Lemma 2. The twist of the fixed point at the origin of the Hamiltonian 
Eq. (flip is given by Eq. (TJj where 

A 3B 2 
t = 2A 

n = ^(8eA-3B 2 ), (13) 

Proof. Upon manipulation of Eq. (0), we can reduce the integrals to calcu- 
late to 

r (h) = ±i h -^l dz 



2n J y 3 
J'"(h) = ±-l(*?-*\<l. 



2tt / V y 5 y 3 , 

To evaluate the twist we need to compute these integrals at h = 0. Since the 
cycles reduce to loops around the origin, this amounts to computing residues 
of the integrands. □ 
A similar calculation for n > 3 in Eq. (|]), shows that the lowest order 
twist To is independent of B. Thus it vanishes only if the bare twist term, 
A = 0. This shows that the tripling normal form is the only one for which 
the twistless torus is generically created at the origin. The reason that the 
tripling is different is that the order of Q(I) is smaller than that of the twist 
term AI 2 . 

Lemma |2| implies that the twist vanishes at the origin when 
3B 2 

e = e = — , (14) 
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which corresponds to the frequency, in the nonrotating frame, of 



1 W 



2 



^0 = 77 + 



3 4A 

The twistless curve moves away from the origin as e moves towards zero, cor- 
responding to approaching the tripling bifurcation. Using these expressions, 
we can compute the rotation number of the twistless curve near e where it 
is created. In the original nonrotating frame, the rotation number is given 
by 

1 A 

Slo{w) =UJ- 77-gj(^ - ^o) 2 + 0{U - UJ ) 3 

where we have expressed it as a function of the rotation number uj of the 
origin. The twistless torus can be seen in Fig. ^, which is a contour plot 
of the rotation number on the bifurcation diagram. Since Q'(h) = at the 
twistless torus, it corresponds to the locus of points where the contours are 
vertical. The twistless torus emerges from the origin at e = 3/4, and collides 
with the saddle-center bifurcation at e = 9/32. In region (ii) the plotted 
curves correspond to the rotation number of the torus in the central triangle; 
in region (iv) they give the rotation number of the torus in the outer region. 
In this way the rotation number can be made continuous. The torus in the 
central triangle corresponding to region (iv) has the negative rotation number 
of the one in the outside region with the same h. This is a result of the sum 

rule £ — = 0, where the ji are the two real cycles for n = 3. 

J 71+72 y ' ' J 

Since the expression for tq is not affected by higher order terms in P or 
Q, we can describe the generic scenario of the tripling bifurcation as follows. 
A sufficient distance from the tripling bifurcation, when e > 3/4, the peri- 
odic orbit has twist, i.e., the rotation number is a monotone function of the 
transverse action near the orbit. Note that this corresponds to u > 1/3 if the 
rotation number increases with action, i.e., A > 0, otherwise it corresponds 
to uj < 1/3. As w moves towards 1/3 or equivalently as e decreases toward 
zero, but before there is any obvious sign of the 1/3 bifurcation, the central 
periodic orbit loses its twist at the parameter value e = 3/4. 

Beyond that point the twist in the island is no longer monotone, and 
there is a twistless torus somewhere in the island. For e < 3/4 a generic 
perturbation of H cannot destroy a diophantine twistless curve with nonzero 
second twist, as ist shown in [ pdlL98|| . The difference to ordinary KAM the- 
ory is that for small perturbations the twistless torus with the same rotation 
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e 

Figure 3: Lines of constant rotation number spaced equidistant with AO = 
0.02 on the bifurcation diagram for the tripling bifurcation. The contours cluster 
around the line of the unstable 3-fold orbit, because there the period T diverges 
logarithmically. 

number exists not necessarily for the same but for slightly change param- 
eter e. At e = 9/32 the twistless torus reaches rotation number 1/3 and 
there is a saddle-center bifurcation creating a pair of period three orbits at 
nonzero radius in the island. The twist inside the central triangle is opposite 
to that outside the period three island chain. Beyond this point the twistless 
torus is replaced by the separatrix with rotation number 1/3. Of course, the 
separatrix is generically a homoclinic tangle. 

Eventually, when u; = l/3or£ = 0, the saddle period three orbit collides 
with the central periodic orbit. Moving past this point, the saddle period 
three orbit re-emerges, and the origin has changed the sign of its twist. Thus 
for e < 0, the rotation number is again a monotone increasing function. 

We conclude by stating this more formally 

Theorem 3. Near a tripling bifurcation of a periodic orbit, a two degree 
of freedom Hamiltonian system has a codimension one twistless bifurcation, 
corresponding to the vanishing of the twist at the periodic orbit. For the 
normal form, Eq. ([U\), this results in the creation of a twistless curve at 
Eo = 3/4, Eq. (\m), that moves away from the origin as the tripling point is 
approached, and eventually collides with the saddle- center bifurcation of the 
period three orbits. 
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3 Tripling bifurcation for area-preserving maps 



In this section we study the normal form for an area-preserving map, /, in 
the neighborhood of an elliptic fixed point with a tripling bifurcation. We 
will show that the twist generically vanishes in the neighborhood of this 
bifurcation. 

Suppose that / has an elliptic fixed point at the origin with rotation 
number u, i.e., with linear multipliers A and A, Eq. (^j). We can put the lin- 
earization of the map in the complex, diagonal normal form (z, z) i— > (Az, Xz) 
by introducing a pair of complex variables (z, z) defined by the linear trans- 
formation 

rp 1 

( x >y) = —{vz-vz) (15) 

where v is the complex eigenvector of Df(0) associated with A. If we choose 
the normalization of the eigenvectors so that v x v = 2i, then the z variables 
are related to action-angle variables by Eq. (||). In the new variables, the 
map has the power series 



A 



3 ?n~3 

n=2 j=0 



+ 0(r + l) (16) 



where the terms are ordered as homogeneous monomials of degree n in z 
and z. The map for z is simply the complex conjugate of Eq. fll~6|) . Note 
that the complex coefficients, not all independent when the map is 

area-preserving. 

The normal form for Eq. flTED can be computed by applying a sequence 
of canonical transformations to eliminate as many terms in the series as 



possible | JSM71jJ . As is well known each of the coefficients can be formally 
eliminated providing the multiplier does not satisfy a "resonance" condition 
of the form 

X j - k ~ 1 = 1 (17) 

In particular for any A, the terms with j — k = 1 cannot be eliminated — 
these we call the "twist" terms, since they give rise to shear in the rotation 
about the fixed point. These terms give rise to the Birkhoff normal form 
Eq. (E3). However, if A is a root of unity, then other resonant terms occur. 
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For example, when A = ~, then we cannot eliminate terms for which j — k — 1 
mod 3. In this case the resonant normal form is 



z 



X [z + a 02 z 2 + a 2 iz 2 z + ai 3 zz 3 + a 40 z A + a 32 z 3 z 2 + ao^z 5 ] + 0(6) 

To study the dynamics in the neighborhood of the tripling bifurcation, 
we assume that the lowest order resonant term, ao 2 , is nonzero. In this case 
we can scale z to eliminate this coefficient. When the map is area-preserving, 
the thirteen real coefficients in the resonant normal form must satisfy four 
conditions through fifth order. After applying these conditions and scaling 
the map to eliminate the coefficient ao 2 , we can reduce the resonant normal 
form to a seven parameter family through 5th order 

z^X[z + z 2 + (l + ia)z 2 z+(f3 + i~f)zz 3 (18) 
+ a i0 z 4 + (5 + ir])z 3 z 2 + a 05 z 5 ] + 0(6) 

The real coefficients a, (3, 7, 77 and the complex coefficient 005 are arbitrary 
and 

a 40 = \[2 + (3 + l (2a- 1 )} , 5 = 2/5 - 



Taking the resonant form Eq. (|T^) as our model, we imagine that the 
coefficients are fixed, and that uj is the bifurcation parameter. When uj 7^ 1/3 
the resonant form can be transformed to Birkhoff normal form. In fact if uj is 
not 1/4, 1/5 or 2/5, we can reduce the map to Birkhoff normal form through 
fifth order. Our goal is to compute the twists To and t% for this map. 

To carry out the transformation we must apply a sequence of four, near 
identity coordinate transformations to eliminate successively all of the quadratic, 
cubic, quartic and quintic terms except for the twist terms. One way to ob- 
tain area-preserving coordinate transformations is to use canonical generating 
functions of the form 

Pq + F {r \P,q) 

where F^' is homogeneous of degree r + 1. This implicitly generates a canon- 
ical transformation through the equations 

p = P + -o r (P,v) , Q = q + - w -{p,q) 
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The lowest order terms of the transformation (q,p) i— > (Q,P) are easily 
obtained explicitly 



P = P- "^-(P. 9) + °( 2r ~ 1). Q = ( l + -g^(P> <?) + 0(2r - 1) , 

(19) 

and the higher order terms can be obtained with more effort, order-by-order. 
Converting the transformation to complex coordinates, with Z = P + iQ, 
yields a transformation of the form 



izV - * + 0(2r - 1) 



where the coefficients by can be explicitly computed in terms of the coeffi- 
cients of F^ r \ 

We wish to compute the nonresonant Birkhoff normal form for Eq. fll8|) 
through fifth order. To do this, we first apply the canonical transformation 
F^ 2 \ with the coefficients chosen to eliminate the quadratic term in Eq. flT8|). 
This is possible whenever u ^ 1/3. This transformation generates terms of 
cubic and higher order in the map. Whenever l/4we can then apply a 
transformation generated by to eliminate all of the cubic terms, except 
for the twist term z 3 z 2 . This transformation will not modify the quartic 
terms in the map, but because of the 0(5) correction terms in Eq. ([19]) , it 
will modify the quintic terms. Finally, whenever u ^ 1/5 or 2/5, we can 
apply transformations F^ and F^ to eliminate all of the quartic and all 
of the quintic terms except for the twist z 3 z 2 . These transformations will 
not modify the coefficient of the twist term, however, and therefore we do 
not need to compute them in order to find the Birkhoff normal form to fifth 
order. This yields the Birkhoff normal form for Eq. ( |i8|) in the form of Eq. (0) 
where 

(3t 2 - 1) 

tttq = a 



t(t 2 - 3) 

vrr, = 4(77 - 7) + 12 ^ ~ \ ( ^ 9 + + vrr ^ ~ \ - (3 ) (20) 
Kl n t(t 2 -3) V2t 2 (t 2 - 3) 2 °t(t 2 - 3) 1 

where t = tan(7ra;). 
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Near e = Owe can reduce these expressions to 
ttto = --— + a + 0(e) 

37T£ 

™ = ^ + 4? T "-h 3+ ^"'- 2T " )+ ° {El) (21) 

The dominant terms in these expressions, though they look quite different, 
are actually equivalent to those obtained for the flow in Lemma |2Q providing 
we set (3 = rj = 7 = 0. Nonzero values would correspond to higher order 
terms in the flow Hamiltonian. 

The function tq(uj) in Eq. (|20|) is shown in Fig. f|. It is monotone increas- 
ing and maps the domain [|, |) U (|, |) one-to-one onto (—00, 00). Moreover, 
the coordinate transformations that we carried out to compute tq are valid 
in this domain. Thus for any a, there is a rotation number ujq for which 
T o{ u) o) = 0- If ^0 is n °t 1/5,1/4 or 2/5 then the transformations leading 
to the expression for T\ are valid. Whenever Ti(uq) 7^ 0, as is generically 
true for our expressions, the Moser twist theorem ||SM71|| implies that there 



are invariant circles in the neighborhood of the elliptic point. Thus we can 
conclude 

Theorem 4. Let f w be area-preserving map with an elliptic fixed point that 
has rotation number u. Suppose that when f is transformed to its resonant 
normal form, Eq. jjT^), the coefficient a is finite (i.e. that a 02 7^ 0). Then 



there is an u £ [|, |) U (|, |) such that / Wo has a twistless bifurcation, i.e. 
where tq(u ) = 0, from Eq. The elliptic fixed point is stable whenever 

ujq $l {1/5, 1/4,2/5}, and Ti(uo) 7^ 0, as is generically the case. 



Note that in Eq. ( |21| ) the dominant contribution to t\(loq), when luq is 
reasonably close to |, has the sign of e. Thus there are two typical scenarios 
for the twistless bifurcation. If uq < |, then typically T\{uq) < 0, and so 
the rotation number Q(J) is locally monotone decreasing as u — > ujq from 
below. The twistless curve is created at a local maximum of Q(J) which 
moves away from the origin as u approaches |. Conversely, if u > |, 

2 The time 2tt map of the flow of Eq. ([ll]), is approximately 

C^x{c + 2mA\C\\+^ C 2 

which is equivalent to Eq. (|l8|), upon identifying a = ^7^, where s = (3irB) 2 /2 is the 
factor that scales the action to normalize 002 in the map. 
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Figure 4: Plot of the twist tq{u) in Eq. (^) with a = 0. 



then typically, ti(uq) > 0, so the rotation number is monotone increasing as 
uj — > cj from above. The twistless curve corresponds to a local minimum of 
fl( J) which moves away from the origin as u moves towards |. On the other 
hand, it is also possible that Ti(lu ) has the opposite sign of e, in which case 
the bifurcation creates a twistless curve as to moves from ujq in the direction 
away from |. 

In the actual dynamics we expect that when the rotation number of the 
twistless curve passes through rational points, there will be reconnection bi- 
furcations [ |HM84| , |HH95| | in which two island chains with the same rotation 



number annihilate each when they collide at the twistless "curve." We ex- 
pect this to happen often for the island chains with rotation numbers in the 
interval [uq, |), because after birth, they will travel outward from the fixed 
point, and if they move rapidly enough will cross the twistless curve. 

In the next sections, we will give examples of the twistless bifurcation for 
the Henon map, and for a quartic map, where we can choose parameters so 
that loq is a low order rational. 

3.1 Tripling bifurcation for the Henon map 

Any quadratic, area-preserving map of the plane can be written in the form 
[Hen69|| 

(x, y) t— > (y — k + x 2 , —x) . 
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It has an elliptic fixed point at x e = —y e = 1 — y/1 + k when — 1 < k < 3. 
The rotation number at the fixed point is 

1 /1 + ^ 1/4 



uj = — arcsm 



7T \ 4 

The normalized, complex eigenvector is v — -^(— A, 1) T , where s = sin(27rc<j). 
Using Eq. (|15D, shifted to the elliptic fixed point, we define the complex 
variable 

z = — -= (x - x e + X(y - y e )) , 

so that the Henon map becomes 
1 



Xz 



Applying the normal form transformation reduces the map to the Birkhoff 
normal form, Eq. (Q), where 

1 (3t 2 -5)(l+t 2 ) 3 
2 6 tt t 4 (3 - 1 2 ) 
1 (1 + t 2 ) 6 (51t 10 - 637t 8 + 2038t 6 - 2706t 4 + 2055t 2 - 705) 



Tl 



2 13 7T t 9 (t 2 - l)(t 2 - 3) 3 

where t = tan(7rco>), as before. Moser previously obtained r in ||Mos94j| P|. 



Note that r vanishes only when t = ■^/5/3, which occurs when u = uJq 
^ arccos(— |)) w 0.2902153116, or k = 9/16. As Moser noted, u is a tran- 
scendental number, and therefore satisfies a Diophantine condition. Accord- 
ing to a theorem of Riissmann, this implies that there are invariant curves 
in the neighborhood of the fixed point so it is stable. This result is obtained 
more easily by noting that 7~i(u;o) ~ —.3366 is nonzero. In this case the 



stability of the fixed point follows from the Moser twist theorem [ |5M71 . 

Both twist functions are shown in Fig. [5] as a function of u. The Birkhoff 
normal form has a curve of zero twist that is born at u and moves outward 
in action as u increases at least up to 0.31375 where T\ vanishes. This 
upper endpoint is artificial, however, since the higher order twists will become 



3 Moser's equivalent expression is apparently not written in standard canonical vari- 
ables, and differs from ours by a scaling factor (of —1/2) 
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Figure 5: Plot of tq (thick line) and t\ (thin line) for the Henon map as a function 
of uj. Here t vanishes only at u « 0.2902153116, and n at u « 0.2308206101, 
0.3137515644, and 0.3944381765. 

increasingly important as n nears zero. The rotation number of the twistless 
curve is approximately 

n (w)»w-^ (22) 

2 Ti 

This function is shown in Fig. |6|. Based on the flow results, we expect that 
the twistless curve in some sense collides with the saddle-center point of the 
period three orbit, k = 1 — this corresponds to uj = u 5 in the figure, so that 
^0(^5) = 1/3- We saw in Sec. |2| that the twistless curve does collide with 
the saddle-center bifurcation for the resonant normal form of the flow, which 
was integrable and for which we obtained an exact expression for O. The 
predictions of the normal form for the map, however, are not of much use 
here because the Henon map exhibits considerable chaos for orbits at these 
radii. If the twistless "curve" continues to the tripling, it will no longer be 
an invariant circle; presumably it will be a cantorus. 

In the interval where the twistless curve exists, u>o < uj < UJ5, there 
should be many unusual bifurcations. For example, whenever the twistless 
curve crosses a rational rotation number, there will typically be a collision 
of two island chains of periodic orbits with this rotation number. The most 
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Figure 6: The rotation number of the twistless circle, £lo(u>) near the tripling 
bifurcation of the Henon map. 



prominent case corresponds to the lowest order rational in this interval, 3/10, 
which we sketch in the phase portraits in Fig. |7|. Similar bifurcations can 
be observed when the twistless curve passes through e.g. 4/13 or 5/16. The 
normal form indicates that the twistless curve has rotation number Qq(u) = 
3/10, when u ~ 0.2995198, using the expansion Eq. (j2~2|) . Translating back to 
the parameter of the Henon map, this corresponds to k m 0.7060175. Indeed 
there is a saddle-center 3/10 bifurcation near k = 0.7063832, which differs 
from our prediction by 0.05%. However, the situation is more complicated 
than the normal form would indicate. We list the bifurcations in Table [l], 
and sketch the corresponding situation in Fig. [8[ The four 3/10 orbits are 
created in two successive 3/10 saddle-center bifurcations at u)\ and uj 2 . In 
the interval uj 2 < uj < there are four 3/10 orbits, two are elliptic and 
two are hyperbolic. The stable and unstable manifolds of the two hyperbolic 
3/10 orbits undergo a reconnection bifurcation | )HH95| | near uj 3 , see Fig. [7]. 
Subsequently the "inner" pair of 3/10 orbits collides with the elliptic fixed 
point, in a decupling bifurcation at 004. = 3/10. The outer pair of 3/10 orbits 
move away from the origin and persist as k — > 00. These last two orbits are 



the orbits that we continue from the "anti-integrable" limit, R5DM98 
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Figure 7: Phase portraits of the Henon map near uj^ where the 3/10 reconnection 
bifurcation occurs, k = 0.70638, 0.70639, 0.70640, 0.70640121, 0.70642, 0.70650 

3.2 Normal form for a quart ic map 

To study codimension two bifurcations, we need a map with two essential 
parameters. Since the quadratic map has only one (the rotation number), we 
turn to a higher order polynomial map which we take to be the composition 
of a rotation and a shear 





| — > R2TTUJ 1 







y + x 2 + cx 3 + dx 4 



where is a counterclockwise rotation by the angle (p. This map has 3 
essential parameters, u, c, and d. According to our general theory, we expect 
that with a choice of one parameter, u, or c we can set r = 0. In fact, with 
the choice of c, we can set u to be any value we choose. Finally, with the 
choice of d we can set t\ = 0, as well. 

Transforming this map into Birkhoff normal form, we obtain Eq. (0), 
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Table 1: Bifurcation Points for the Henon map 



Bifurcation 


label 


UJ 


A 


twistless 




0.2902153 


9 
16 


sn(3/10) 


LOi 


0.2995432 


0.7063832 


sn(3/10) 


UJ 2 


0.2995438 


0.7063926 


reconnection 


UJ 3 


0.2995444 


0.70640121 


decupling 


UJ4 


3 
10 


0.7135255 


sn(l/3) 


U) 5 


0.3179717 


1 


tripling 


UJ 6 


1 
3 


5 
1 



providing the nonresonance conditions are satisfied as usual. The twists are 
7ir 



3c 5 - 3t 2 



t(t 2 - 3) 

-705 + 2055t 2 - 2706t 4 + 2038t 6 - 637t 8 + 51t 10 



3c 



+ 



256t(t 2 - 1) 
3d 7-5t 2 



512t 3 (t 2 - 3) 3 (t 2 - 1) 

-225 + 556t 2 - 534^ + 204t 6 -m 8 c 

t(t 2 -3) 2 2 V ; 



16 t(t 2 - 3) 

For the cubic map with d = 0we can make To = at any ujq 7^ 1/3 with an 
appropriate choice of c. For example, we can set u = 1/5 or 2/5, and then 
choose 



6 



10±2V5. 



When r = at a fifth order resonance, Th. |5| in the Appendix implies that 
the fixed point is unstable at the bifurcation point. The corresponding map 
with u ~ 1/5 is shown in Fig. |9|. The instability is difficult to see numerically 
because the size of the nonlinear terms is so small near the fixed point, though 
in the figure it is clear that the unstable motion does approach the origin. 
It is interesting to note that the instability is much more pronounced for 
uj = 1/5 than for u = 2/5. 
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Figure 8: Sketch of the 3/10 and 1/3 bifurcations for the Henon map as a function 
of u. Representative phase portraits are also shown. See Table |l| for bifurcation 
values. 

The second twist T\ can be used in two different ways. For the cubic map 
with d = we can show that if tq = then always T\ ^ 0, so that the only 
exceptions to stability can be the third through sixth order resonances. The 
value of c that makes tq = at t = to and the corresponding value of T\ are 

5 - 3t 2 _ 105 - 305t| + 353t^ - 167^ + 30t§ 

C "3to(tg-3)' 7fTl - 48lf(tl - 3)3(tg - 1) • 

It is easy to see that the polynomial in the numerator of T\ is positive, so 
that T\ is always nonzero if r vanishes for the cubic map. 

For the quartic map with d ^ we can make both twists vanish at 
arbitrary Uq, in particular at seventh order resonance. In general 

5-3t 2 105 - 305tg + 353^ - 167^ + 30tg 

C "3t (tg-3)' 9*g(*g - 3)=(^ - l)(5t= - 7) ' 
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and with t = tan(vr/7) we obtain c « 1.0763012774 and d « .7144291292. 
By Th. |5| the fixed point is unstable for these parameter values. Note that 
if we let c and (i be the above functions of uj we obtain the very special one 
parameter family of quartic maps for which the two first twists vanish for 
every frequency. 




Figure 9: Phase portrait of the cubic map (d=0) with ro = Oatu; = l/5 — 
e, 1/5, 1/5 + e, e = 0.0004. At uj = 1/5 the fixed point is shown to be unstable. 



4 Conclusion 

By analyzing the normal form for Hamiltonian flows and area-preserving 
maps in the neighborhood of the tripling bifurcation we have shown that 
there exists a bifurcation creating a twistless torus. Thus in the presence of 
a tripling bifurcation, the non-degeneracy condition of the KAM theorem is 
violated on a one parameter family of tori in phase space in the neighbor- 
hood of the bifurcation. Moreover, in the nonintegrable case, reconnection 
bifurcations occur for orbits that collide at the twistless torus, such as the 
3/10 orbits of the Henon map in Fig. |8]. 

Using the notation of Fig. §, we say that a twistless bifurcation occurs 
at uj = uj . The pair of period three orbits is created in a saddle-center 
bifurcation at uj$, and the tripling occurs at ujq = 1/3. For the case of 
the resonant Hamiltonian normal form, Eq. (|TJ]) implies that the relative 
positions of these bifurcations have a simple ratio 

uj - 1/3 _ 8 
uj 5 — 1/3 3 
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Since the Hamiltonian resonant normal form is integrable, we can follow the 
twistless torus; it exists when the rotation number of the central periodic or- 
bit is in the interval ujq to uj^. Moreover, the rotation number of the twistless 
curve itself, Qq, takes all of the values in the interval from ujq to 1/3. For 
every rotation number in this interval there will generically be a reconnection 
bifurcation for the full two degree of freedom Hamiltonian. Hence this type 
of bifurcation occurs generically in nonintegrable Hamiltonian systems. 

Similarly any area-preserving map that has a tripling bifurcation of an el- 
liptic fixed point and that satisfies the nondegeneracy condition that the third 
order resonant coefficient, 002 does not vanish, has a twistless bifurcation. 
The dependence of the twist on rotation number is more complicated for this 
case than for the Hamiltonian, nevertheless, the rotation number of the fixed 
point unfolds the bifurcation. Moreover, there is always an w 6 [1/6,1/2) 
for which the resonant normal form has a twistless bifurcation. 

As an example, the Henon map has a twistless bifurcation at k — 9/16, 
and the twistless curve moves outward as k increases, causing reconnec- 
tion bifurcations. We examined in particular the bifurcations of the 3/10 
orbits. The vanishing of twist (s), as we showed for a cubic/quartic area- 
preserving map, can also lead to the instability of the fixed point at quintu- 
pling/ septupling bifurcations. 

The whole scenario is similar for the quadrupling bifurcation, when the 
resonant term dominates the twist term in magnitude (see the Appendix) . In 
this case a pair of period 4 orbits are created in a saddle-center bifurcation, 
and a twistless torus is created at a cusp in the bifurcation diagram similar 
to that in Fig. [l]. However, in this case the twistless curve does not collide 
with the fixed point, and a simple prediction about what rotation numbers 
must pass through the twistless curve is not possible. 

It would be interesting to generalize these results to higher dimensional 
systems. 

5 Appendix: 

Instability of resonant maps without twist 

In this section we give sufficient conditions for which a map of the plane with 
eigenvalues that are n th roots of unity has an unstable fixed point. Here the 
condition of area preservation is not necessary. Examples of this behavior are 
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well known [ |SM71|| , however, we have not found an explicit proof of instability 
for the case where the resonant term has the lowest possible order. 

To study this problem, we put the map in the normal form Eq. (^). We 
will show that the fixed point is generically unstable whenever all of the twist 
terms vanish that are of lower order than 0(n— 1), the order of the resonant 
term. When n is even, there is a twist term of the same order as the resonant 
term. In this case the 0(n— 1) twist term need not vanish, but its magnitude 
must be dominated by the resonant term. 

Theorem 5. Let f be a C n map of the plane with a fixed point at the origin. 
Suppose that the multipliers of the fixed point are n th roots of unity, and that 
when f is put in normal form, the low order twist coefficients vanish: 

n 

= > 1 < 3 < j ' 
but the lowest order resonant term is nonzero 

O0,n-1 7^ • 

If n is even, assume in addition that 

\a>n/2,n/2-l\ < \a>0,n-l\ ■ 

Then the origin is unstable. 

Note that when n = 3, the assumed form is generic, but for n = 4, the 
first twist, a 2 i, is assumed to be small, and when n = 5 it is assumed to 
vanish. 

The instability at the (generic) tripling bifurcation is well known. Our 
calculation shows that a similar phenomenon occurs for the quintupling bi- 
furcation: the fixed point becomes unstable if the first twist vanishes. As we 



saw in Sec. |3.2| this is a generic bifurcation in a two parameter area-preserving 
family, one parameter to adjust the frequency and the second one to make 
the first twist vanish. 

Proof. When n is odd the assumed normal form is 

z' = \{z + az n - l )+0(n) . 



We adapt the proof from Siegel and Moser ||SM71|| . Begin by scaling z using 
the transformation z —>■ zs where s determined by s n_1 /s = (net) -1 . Then, 
using the fact that A™ = 1, the n th power of z' becomes 

( z ') n = z n + z n - x z n ~ x + 0(2n - 1) . 
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Introducing new coordinates Z by 
z n = Z = X + iY = i?expi$ , 

gives 

z' = z + R 2 - 2/n (1 + n(Z, Z)R 1/n ) , 

where the higher order terms have been represented by a factor 77, which, 
since / is C n , is bounded as R — > 0. The purpose of the scaling to eliminate 
a is that the real part of this map now becomes 

X' = X + R 2 ~ 2/n (1 + ${r])R 1/n ) . 

Since rj is bounded, there is an R Q such that for all R < R$, we have 
\$t(?l)R 1/n \ < 1/2. Thus whenever R < R , X' > X + \R 2 ~ 2 l n , and so 
X is monotone increasing in each step of the map. It is easy to see that this 
implies that the fixed point is unstable. 

The case of even n is slightly harder because there can be a twist term of 
the same order: 

z' = A (z + az n ~ x + (3z\z\ (n - 2) + 0{2n - 1)) . 

If f3 happens to be zero the above proof works. When f3 is nonzero the 
stability of the fixed point is determined by the relative size of the resonant 
term and the twist term. Again we scale away a and the n th power of the 
new map is 

( Z >Y = z n + (zzy- 1 + (3z n \z\( n -V + 0(2n - 1), 
where (3 — /3/\a\. In the new coordinates this becomes 

Z' = Z + R 2 - 2/n (l + exp(i$) + r]R 1/n ), 
and the map for the real part is 

X' = X + R 2 - 2 ' n (\ + K(/3exp(i$)) + ^(ri)R 1 ^ . 

So if 1 + 3fJ(/9exp(i$)) > for all $ the argument of the previous case 
applies. Therefore the instability criterion is |/3|/|o;| < 1: if the "twist" term 
is smaller than the resonant term the mapping is unstable. The well known 
example of this type is the quadrupling bifurcation. Our calculation shows 
that a similar codimension two bifurcation occurs for the area preserving 
sextupling bifurcation. □ 
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